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ABSTRACT 

We summarize our calculation of the total cross section for 
top quark production at hadron colliders within the context of 
perturbative quantum chromodynamics, including resummation 
of the effects of initial-state soft gluon radiation to all orders in 
the strong coupling strength. 

I. INTRODUCTION AND MOTIVATION 

In hadron interactions at collider energies, tt pair production 
proceeds through partonic hard-scattering processes involving 
initial-state light quarks q and gluons <?. In lowest-order pertur- 
bative quantum chromodynamics (QCD), 0{a s ), the two par- 
tonic subprocesses are q + q t + i and g + g — > t+i. Calcula- 
tions of the cross section through next-to-leading order, O(o%), 
involve gluonic radiative corrections to these lowest-order sub- 
processes as well as contributions from the q + g initial state 
A complete fixed-order calculation at order 0(a™), n > 4 does 
not exist. 

The physical cross section for each production channel is ob- 
tained through the convolution 



aij(S, m) 



4m 2 



dq$ i j{r),n)&ij{r),m,n). (1) 



The square of the total hadronic center-of-mass energy is S, 
the square of the partonic center-of-mass energy is s, m de- 
notes the top mass, fj, is the usual factorization and renormal- 
ization scale, and is the parton flux. The variable 

V = — 1 measures the distance from the partonic thresh- 
old. The indices ij £ {qq, gg} denote the initial parton channel. 
The partonic cross section <xy(?7,m, fi) is obtained commonly 
from fixed-order QCD calculations or, as described here, 
from calculations that go beyond fixed-order perturbation the- 
ory through the inclusion of gluon resummation [Q, ||, [I]] to all 
orders in the strong coupling strength a s . We use the notation 
a = a(fj, = m) = a s (m) / tt. The total physical cross section is 
obtained after incoherent addition of the contributions from the 
the qq and gg production channels. 

Comparison of the partonic cross section at next-to-leading 
order with its lowest-order value reveals that the ratio becomes 
very large in the near-threshold region. Indeed, as rj — > 0, the 
"Jf-factor" at the partonic level K{rj) grows in proportion to 
o In 2 (77). The very large mass of the top quark notwithstanding, 
the large ratio K(rj) makes it evident that the next-to-leading 
order result does not necessarily provide a reliable quantitative 
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prediction of the top quark production cross section at the en- 
ergy of the Tevatron collider. The large ratio casts doubt on the 
reliability of simple fixed-order perturbation theory for physi- 
cal processes for which the near-threshold region in the suben- 
ergy variable contributes significantly to the physical cross sec- 
tion. Top quark production at the Fermilab Tevatron is one such 
process, because the top mass is relatively large compared to 
the energy available. Other examples include the production of 
hadronic jets that carry large values of transverse momentum 
and the production of pairs of supersymmetric particles with 
large mass. To obtain more reliable theoretical estimates of the 
cross section in perturbative QCD, it is important first to iden- 
tify and isolate the terms that provide the large next-to-leading 
order enhancement and then to resum these effects to all orders 
in the strong coupling strength. 

II. GLUON RADIATION AND 
RESUMMATION 

The origin of the large threshold enhancement may be traced 
to initial-state gluonic radiative corrections to the lowest-order 
channels. We remark that we are calculating the inclusive total 
cross section for the production of a top quark-antiquark pair, 
i.e., the total cross section for t + i + anything. The partonic 
subenergy threshold in question is the threshold for t + i+ any 
number of gluons. This coincides with the threshold in the in- 
variant mass of the t+i system for the lowest order subprocesses 
only. 

For i + j — > t + i + g, we define the variable z through the 
invariant (1 — z) = f * , where k and p t are the four- vector 
momenta of the gluon and top quark. In the limit that z — > 1, 
the radiated gluon carries zero momentum. After cancellation 
of soft singularities and factorization of collinear singularities 
in 0(atg), there is a left-over integrable large logarithmic con- 
tribution to the partonic cross section associated with initial- 
state gluon radiation. This contribution is often expressed in 
terms of "plus" distributions. In O(o%), it is proportional to 
a 3 ln 2 (l — z). When integrated over the near-threshold region 
1 > z > 0, it provides an excellent approximation to the full 
next-to-leading order physical cross section as a function of the 
top mass. At m = 175 GeV, the ratio of the next-to-leading 
order to the leading order physical cross sections in the lead- 
ing logarithmic approximation is cr^ +1 /cr§i =1.22. This ra- 
tio shows that the near-threshold logarithm builds up cross sec- 
tion in a worrisome fashion. It suggests that perturbation the- 
ory is not converging to a stable prediction of the cross sec- 
tion. The goal of gluon resummation is to sum the series in 
a n+2 i n 2n (i _ z j t a jj orders in a in order to obtain a more 
defensible prediction. 

Different methods of resummation differ in theoretically and 



phenomenologically important respects. Formally, if not explic- 
itly in some approaches, an integral over the radiated gluon mo- 
mentum z must be done over regions in which z — * 1. There- 
fore, one significant distinction among methods has to do with 
how the inevitable "non-perturbative" region is handled. In 
the approach of Laenen, Smith, and van Neerven (LSvN) [[|], 
an undetermined infrared cutoff (IRC) p, is introduced, with 
A-QCD < Po < fn. The presence of an extra scale spoils the 
renormalization group properties of the overall expression. The 
unfortunate dependence of the resummed cross section on this 
undetermined cutoff is important numerically since it appears in 
an exponent [|p. It is difficult to evaluate theoretical uncertain- 
ties in a method that requires an undetermined infrared cutoff. 



III. PERTURBATIVE RES UMM ATION 

The method of resummation we employ [Q| is based on a per- 
turbative truncation of principal-value resummation [^J. This 
approach has an important technical advantage in that it does 
not depend on arbitrary infrared cutoffs. Because extra scales 
are absent, the method permits an evaluation of its perturba- 
tive regime of applicability, i.e., the region of the gluon radia- 
tion phase space where perturbation theory should be valid. We 
work in the MS factorization scheme. 

Factorization and evolution lead directly to exponentiation of 
the set of large threshold logarithms in moment (n) space in 
terms of an exponent E PV : 



E PV (n,m 2 ) 



n-l 



dX 

i-C J ~ 

(i-C) 2 



g [a (Am 2 )] . 



(2) 

The function g(a) is calculable perturbatively. All large soft- 
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Figure 1 : Differential cross section da / dij in the MS scheme for 
the qq channel: Born (dotted), next-to-leading order (dashed) 
and resummed (solid). 



gluon threshold contributions are included through the two-loop 
running of a. The integral in the complex plane runs along 
a contour P with endpoints and 1 that is symmetric under 
reflections across the real axis. 
The function E pv is finite, and lim„_»oo E PV (n, m 2 ) 
x . Therefore, the corresponding partonic cross section is fi- 
nite as z — > 1 (n — > +oo). The function E pv includes both 
perturbative and non-perturbative content. The non-perturbative 
content is not a prediction of perturbative QCD. We choose to 
use the exponent only in the interval in moment space in which 
the perturbative content dominates. We use the attractive finite- 
ness of Eq. (Q) to derive a perturbative asymptotic representation 
of E(x, a(mj) that is valid in the moment-space interval 



1 < x = Inn < t — 



1 



2a& 2 



(3) 



The coefficient 62 = (UCa — 2n/)/12; the number of flavors 
nf = 5; C q q = Cp = 4/3; and C gg = Ca = 3. 
The perturbative asymptotic representation is 



N(t)+1 p+1 



Eij (x, a) ~ Eij (x, a, N(t)) = 2C tj ^ a " ' 



p=i 3=0 



Here 



(4) 



(5) 



and T(l + z) = X^fcLo c k zk > where T is the Euler gamma 
function. The number of perturbative terms N(t) in Eq. (Q) 
is obtained [[}J by optimizing the asymptotic approximation 



E(x,a) - E(x,a,N(t)) 



minimum. Optimization works 



perfectly, with N(t) — 6 at m — 175 GeV. As long as n is in 
the interval of Eq. (^|), all the members of the family in n are 
optimized at the same N(t), showing that the optimum number 
of perturbative terms is a function of t, i.e., of rn only. 

Because of the range of validity in Eq. ([}]), terms in the ex- 
ponent of the form a k ln fc n are of order unity, and terms with 
fewer powers of logarithms, a k ln fe ~ m n, are negligible. Re- 
summation is completed in a finite number of steps. Upon using 
the running of the coupling strength a up to two loops only, all 
monomials of the form a k ln fe+1 n, a k ln fc n are produced in 
the exponent of Eq. (Q). We discard monomials a k ln fc n in the 
exponent because of the restricted leading-logarithm universal- 
ity between tt production and massive lepton-pair production, 
the Drell-Yan process. 

The exponent we use is the truncation 



E ij (x,a,N)=2C i 



N(t) + 1 

E 

P =i 



a"s x p+ \ 



(6) 



with the coefficients s r 



s p+x , p = b^P/pip + 1). The 
number of perturbative terms N(t) is a function of only the top 
quark mass m. This expression contains no factorially-growing 
(renormalon) terms. It is valuable to stress that we can derive 
the perturbative expressions, Eqs. (Q), (Q), and (||), without the 
principal-value prescription, although with less certitude [BR. 



After inversion of the Mellin transform from moment space 
to the physically relevant momentum space, the resummed par- 
tonic cross sections, including all large threshold corrections, 
can be written 



- R:pert 



dze^^^'^a'iMm^). (7) 



The leading large threshold corrections are contained in the ex- 
ponent Eij(x,a), a calculable polynomial in x. The deriva- 



tive <7y (77, to, z) — d(alj' (77, to, z))/ dz, and is the lowest- 
order 0(a 2 s ) partonic cross section expressed in terms of in- 
elastic kinematic variables. The upper limit of integration, 
Zmax < 1, is set by the boundary between the perturbative and 
non-perturbative regimes, well specified within the context of 
the calculation, and z m i n is fixed by kinematics. 

Perturbative resummation probes the threshold down to 77 > 
770 = (1 — Zmax)/^- Below this value, perturbation theory, re- 
summed or otherwise, is not to be trusted. For to = 175 GeV, 
we determine that the perturbative regime is restricted to val- 
ues of the subenergy greater than 1.22 GeV above the threshold 
(2m) in the qq channel and 8.64 GeV above threshold in the gg 
channel. The difference reflects the larger color factor in the gg 
case. The value 1.22 GeV is comparable to the decay width of 
the top quark. 

IV. PHYSICAL CROSS SECTION 

Other than the top mass, the only undetermined scales are the 
QCD factorization and renormalization scales. We adopt a com- 
mon value ft for both, and we vary this scale over the interval 
/i/m <E {0.5, 2} in order to evaluate the theoretical uncertainty 
of the numerical predictions. We use the CTEQ3M parton den- 
sities [Bp. A quantity of phenomenological interest is the dif- 
ferential cross section d<T *i( s ^™ ^ m jt s integral over 77 is the 
total cross section. In Fig. 1 we plot this distribution for the qq 
channel at to = 175 GeV, y/S = 1.8 TeV, and fi = to. The 
full range of 77 extends to 25, but we display the behavior only 
in the near-threshold region where resummation is important. 
We observe that, at the energy of the Tevatron, resummation is 
quite significant for the qq channel. A similar figure for the gg 
channel may be found in our publications [^J. 

In Fig. 2, we show our total cross section for /j/5-production 
as a function of top mass in pp collisions at y/~S =1.8 TeV. 
The central value is obtained with the choice fi/m = 1, and 
the lower and upper limits are the maximum and minimum of 
the cross section in the range /j/m 6 {0.5,2}. At m = 175 
GeV, the full width of the uncertainty band is about 10% . We 
consider that the variation of the cross section over the range 
n/rn £ {0.5, 2} provides a good overall estimate of uncertainty. 
For comparison, we note that over the same range of /j, the 
strong coupling strength a varies by ±10% at to = 175 GeV. In 
estimating uncertainties, we do not consider explicit variations 
of the non-perturbative cutoff, expressed through z max . This 
is justified because, for a fixed to and \i, z max is obtained by 
enforcing dominance of the universal leading logarithmic terms 
over the subleading ones. Therefore, z max is derived and is 
not a source of uncertainty. At fixed to, the cutoff necessarily 



varies as /j and thus a vary. Using a different choice of parton 
densities [PJj, we find a 4% difference in the central value of our 
prediction [Bp at to = 175 GeV. A comparison of the predic- 
tions |Q| in the MS and DIS factorization schemes also shows a 
modest difference at the level of 4%. 

Our calculation is in agreement with the data [p|. We find 
cr**(m = 175 GeV,y/S = 1.8 TeV) = 5.52±g;^ pb. This 
cross section is larger than the next-to-leading order value by 
about 9%. 

The top quark cross section increases quickly with the en- 
ergy of the pp collider. We provide predictions in Fig. 3 for an 
upgraded Tevatron operating at yfS = 2 TeV. We determine 
cr"(m = 175 GeV, y/S = 2 TeV) = 7.56±g;^ pb. The 
central value rises to 22.4 pb at y/S = 3 TeV and 46 pb at 
y/S = 4 TeV. 

Extending our calculation to much larger values of m at 
y/S = 1.8 TeV, we find that resummation in the principal qq 
channel produces enhancements over the next-to-leading order 
cross section of 21%, 26%, and 34%, respectively, for m = 
500, 600, and 700 GeV. The reason for the increase of the en- 
hancements with mass at fixed energy is that the threshold re- 
gion becomes increasingly dominant. Since the qq channel also 
dominates in the production of hadronic jets at very large values 
of transverse momenta, we suggest that on the order of 25% of 
the excess cross section reported by the CDF collaboration ^ 
may well be accounted for by resummation. 

Turning to pp scattering at the energies of the Large Hadron 
Collider (LHC) at CERN, we note a few significant differences 
from pp scattering at the energy of the Tevatron. The dom- 
inance of the qq production channel is replaced by gg domi- 
nance at the LHC. Owing to the much larger value of y/S, the 
near-threshold region in the subenergy variable is relatively less 
important, reducing the significance of initial-state soft gluon 
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Figure 2: Inclusive total cross section for top quark production 
in pp collisions at y/S =1.8 TeV. The dashed curves show the 
upper and lower limits while the solid curve is our central pre- 
diction. CDF and DO data are shown. 



radiation. Lastly, physics in the region of large subenergy, 
where straightforward next-to-leading order QCD is also in- 
adequate, becomes significant for ti production at LHC ener- 
gies. Using the approach described in this paper, we estimate 
a tl (m = 175 GeV, \/S = 14 TeV) = 760 pb. 

V. OTHER METHODS OF RESUMMATION 

Two other groups have published calculations of the to- 
tal cross section at m = 175 GeV and yfs — 1.8 TeV: 
o-"(LSvN [§) = 4.95±°;™ pb; and ct"(CMNT [§) = 4.75±g;|l 
pb. From the purely numerical point of view, all agree within 
their estimates of theoretical uncertainty. However, the resum- 
mation methods differ as do the methods for estimating uncer- 
tainties. Both the central value and the band of uncertainty of 
the LSvN predictions are sensitive to their arbitrary infrared cut- 
offs. To estimate theoretical uncertainty, we use the standard /i- 
variation, whereas LSvN obtain theirs primarily from variations 
of their cutoffs. 

The group of Catani, Mangano, Nason, and Trentadue 
(CMNT) [Q] calculate a central value of the resummed cross 
section (also with /j,/m — 1) that is less than 1% above the ex- 
act next-to-leading order value. There are similarities and dif- 
ferences between our approach and the method of CMNT. We 
use the same universal leading-logarithm expression in moment 
space, but differences occur after the transformation to momen- 
tum space. The differences can be stated more explicitly if we 
examine the perturbative expansion of the resummed hard ker- 
nel H*j(z, a). If, instead of restricting the resummation to the 
universal leading logarithms only, we were to use the full con- 
tent of Hfj (z, a), we would arrive at an analytic expression that 
is equivalent to the numerical inversion of CMNT, 



1 + 2aCi 



ln 2 (l-z) + 2 7 , B m(l 



0(a 2 ). (8) 



In terms of this expansion, in our work we retain only the lead- 
ing term ln 2 (l — z) at order a, but CMNT retain both this term 
and the subleading term 2^e ln(l — z). Indeed, if the sublead- 
ing term 27sln(l — z) is discarded in Eq. (|8|), the residuals 
<5 y /er^ LO defined by CMNT [§] increase from 0.18% to 1.3% 
in the qq production channel and from 5.4% to 20.2% in the gg 
channel. After addition of the two channels, the total residual 
5/a NLO grows from the negligible value of about 0.8% cited 
by CMNT to the value 3.5%. While still smaller than the in- 
crease of about 9% that we obtain, the increase of 3.5% vs. 
0.8% shows the substantial influence of the subleading loga- 
rithmic terms retained by CMNT. 

We judge that it is not appropriate to keep the subleading term 
for several reasons: it is not universal; it is not the same as the 
subleading term in the exact 0(a 3 ) calculation; and it can be 
changed arbitrarily if one elects to keep non-leading terms in 
moment space. The subleading term is negative, and it is nu- 
merically very significant when integrated throughout the phase 
space. In the qq channel at m = 175 GeV and y/S = 1.8 
TeV, its inclusion eliminates more than half of the contribu- 
tion from the leading term. In our view, the presence of nu- 
merically significant subleading contributions begs the question 



of consistency. The influence of subleading terms is amplified 
at higher orders where additional subleading structures occur 
in the CMNT approach with significant numerical coefficients 
proportional to tt 2 , £(3), and so forth. We will present a more 
detailed discussion of these points elsewhere. 

Our theoretical analysis and the stability of our cross section 
under variation of the hard scale /i provide confidence that our 
perturbative resummation yields an accurate calculation of the 
inclusive top quark cross section at Tevatron energies and ex- 
hausts present understanding of the perturbative content of the 
theory. 
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Figure 3: Inclusive total cross section for top quark production 
in pp collisions at \/S = 2.0 TeV. 



